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Abstract 

P 

We studied firstly the ground state properties of C-isotopes using a deformed 
Hartree-Fock (HF)+ BCS model with Skyrme interactions. Shallow deformation 
minima are found in several neutron— rich C-isotopes. It is shown also that the 
deformation minima appear in both the oblate and the prolate sides in 17 C and 
19 C having almost the same binding energies. Secondly, we carried out shell model 
calculations to study electromagnetic moments and electric dipole transitions of the 
C-isotopes. We point out the clear configuration dependence of the quadrupole 
and magnetic moments in the odd C-isotopes, which will be useful to find out the 
deformations and the spin-parities of the ground states of these nuclei. We studied 
electric dipole states of C-isotopes focusing on the interplay between low energy 
Pigmy strength and giant dipole resonances. Reasonable agreement is obtained 
with available experimental data for the photoreaction cross sections both in the 
low energy region below Hco=14 MeV and in the high energy giant resonance region 
(14 MeV < Tilu <30 MeV). The calculated transition strength below dipole giant 
resonance {Uui <14 MeV) in heavier C-isotopes than 15 C is found to exhaust about 
12 ~ 16% of the classical Thomas- Reiche-Kuhn sum rule value and 50 ~ 80% of the 
cluster sum rule value. 



PACS numbers/21.60.Cs, 24.30. Cz, 25.20.-x 
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1 INTRODUCTION 



The structure of nuclei far from the (3— stability lines is often very different from that 
of stable nuclei due to largely extended wave functions as well as the large asymmetry 
between neutron and proton mean fields. Because of these unique features of the mean 
field, a naive extrapolation of the (3— stable nuclei will fail to predict the structure of these 
exotic nuclei. A typical example is the new shell structure at the neutron number N =16 
in O isotopes JTJ. The structure of dipole excitations in neutron-rich O isotopes is also 
found out very different from that of stable nuclei, especially in the low energy region 
below GDR§@, |. 

A study of heavy C-isotopes is another current issue, where these exotic structures may 
be expected. In this paper, we study the ground state properties of C isotopes performing 
deformed Hartree-Fock (HF) +BCS calculations and also shell model calculations. The 
energy surfaces of C isotopes are studied as a function of the quadrupole deformation 
parameter (3% in order to find out the deformation minimum for different single-particle 
configurations. Special emphasis will be put on the magnetic and quadrupole moments 
(Q-moments) for odd C isotopes which will manifest their exotic structure. Electric dipole 
excitations of C isotopes are investigated by a large scale shell model calculations focusing 
on the interplay between low energy Pigmy strength and giant dipole resonance (GDR). 
We try to find out the deformation effect on the dipole strength distributions, which will 
increase the width of GDR. 

The paper is organized as follows. In section 2, we present results of the deformed 
Skyrme HF +BCS calculations. The magnetic and Q- moments are discussed in section 
3. The Pigmy and GDR are shown in section 4. A summary and conclusions are given in 
section 5. 

2 DEFORMED SKYRME HF CALCULATIONS 

In this section, we investigate the neutron number dependence of deformation properties 
along the chain of C-isotopes. For this purpose, we perform deformed HF+BCS calcula- 
tions with Skyrme interactions SkI4 0, SIII and SkM* 0, using the computer code 
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SKY AX ||. The axial symmetry is assumed for the HF deformed potential. The pair- 
ing interaction is treated in the BCS approximation and is taken to be a volume-type 
5-interaction 



where the pairing strength Vq is taken to be —323 MeV-fm 3 for neutron and —310 MeV-fm 3 
for proton || . These values are determined so as to fit the experimental pairing gaps for 
several isotope as well as isotone chains of semi-magic nuclei, see Ref. for details. The 
HF+BCS calculations are also carried out with a density dependent pairing interaction 



where p(r) is the HF density at r = (ri + r 2 )/2 and po is chosen to be 0.16 fm -3 . The 
pairing strength V ' is taken to be —999 MeV-fm 3 for neutrons and —1146 MeV-fm 3 for 
protons . A smooth energy cut-off is employed in the BCS calculations || . 

Figure 1 shows the binding energy surfaces for even-mass C isotopes as a function 
of the quadrupole deformation parameter fa obtained with the SkI4 interaction together 
with the two pairing interactions. In general, both of the two pairing interactions give a 
similar energy surface which is flat in a wide range of the deformation parameter fa. The 
energy minima are tabulated in Table 1. For the volume pairing, the energy surface of 
12 C is rather flat between —0.3 < fa < 0.3. The energy minimum becomes apparent for 
14 C. For heavier C isotopes 16 C and 18 C, two shallow minima appear both in the prolate 
and oblate sides. In 18 C, the ground state has the largest deformation at fa=0.38, while 
the local minimum appears at the oblate side at fa ~ —0.3. In the case of the surface 
pairing (eq. (Q)), clear minima are not seen in the energy surface except for 14 C, although 
there are large flat plateaus between —0.3 < fa < 0.3 as in the case of the volume pairing 
shown in the upper panel of Fig. 1 . 

The surface pairing tends to yield a larger pairing gap for C isotopes than the volume 
pairing (see Table 1). In order to asses whether the slight difference of the energy surfaces 
obtained with the two pairing interaction is due to the different form of interaction or not, 
we repeat the same calculations for the surface pairing interaction but by reducing the 
strength V ' by a half. Even smaller pairing strength is adopted in the study of O isotopes 



V{r 1 ,r 2 ) = V 6(r 1 -r 2 ) 



(1) 
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recently in ref. [fTOfl . Figure 2 shows the energy surfaces thus obtained. As we see, clear 
minima now appear for 14 C and 16 C, suggesting that the energy surface is sensitive to the 
strength of the pairing interaction. 

The results with two different Skyrme interactions SIII and SkM* together with the 
surface pairing are shown in Fig. 3. Although the strengths of the pairing interaction may 
be different for each set of Skyrme interactions, we use the same values as those used with 
the SkM set. The energy surface for the C-isotopes show no clear minima, but rather flat 
in the deformation region —0.3 < /3 2 < 0.3. These features are similar as those of SkM 
interaction. 

The one quasi-particle state energies obtained with the SkM interaction are shown in 
Figs. 4, 5, 6 and 7 for 13 C, 15 C, 17 C and 19 C, respectively. The odd nuclei 13 > 15 > 17 > 19 C are 
treated as one quasi-particle state on top of the BCS ground state of neighboring even 
nuclei. The Pauli blocking effect of the valence particle is not taken into account in the 
present calculations. For 13 C, the l/2~ state shows the deepest spherical minimum at 
/3 2 = 0.0, while the two minima are seen in the l/2 + state at the oblate and the prolate 
deformations. The 5/2 + state has one minimum at around /3 2 = —0.36. The two pairing 
interactions ([!]) and (fj) give essentially the same results for 13 C. These energy minima for 
the different configurations are expected from the Nilsson diagram of deformed harmonic 
oscillator potential |TT[. The results of 15 C are shown in Fig. 5. The l/2 + state is the 
lowest at the prolate deformation with /3 2 ~0.2 with the volume pairing and /? 2 ~0.12 
with the surface pairing. There is also a local minimum for the l/2 + configuration in 
the oblate side. The minimum of 5/2 + state appears at the oblate side and the energy 
is 2MeV above that of l/2 + state. The three configurations l/2 + , 3/2 + and 5/2 + are 
competing in 17 C as shown in Fig. 6. In the case of the volume pairing, the 3/2 + state 
has the lowest minimum at the prolate deformation /3 2 ~0.4 and another local minimum 
is also found at the oblate side with /3 2 ~ —0.2. The energy surface of the l/2 + state 
is similar to that of 3/2 + state although the oblate minimum is lower than the prolate 
minimum in the case of the l/2 + state. The 5/2 + state has a minimum at the oblate 
deformation. The minima of the three configurations at the oblate side show very similar 
/3 2 values having almost the same binding energies. Among all minima, the lowest one is 
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found at the prolate deformation at fa ~0.4 for the 3/2 + state. The binding energies of 
the three states are also very close in the shell model calculations of 17 C, as we discuss in 
the next section. The results of 19 C are shown in Fig. 7. The competition of the three 
configurations l/2 + , 3/2 + and 5/2 + is apparent as in 17 C. The one quasi-particle states 
show the lowest minimum for the 3/2 + state at the oblate deformation fa ~ —0.3, while 
the l/2 + state show the minimum at the prolate deformation fa ~0.32. The 5/2 + state 
does not show any clear minimum as the one quasi-particle state. 

In general, the prolate deformation occurs at the beginning of the shell while the 
oblate deformation occurs at the end of the shell. Experimentally, 12 C and 13 C are known 



to have oblate defer mat ion [ 12 , p~3f| , and 14 C with iV=8 becomes almost spherical at the 
neutron shell closure. In the present calculations, the energy surface of 12 C is rather flat 
between -0.3 < fa < 0.3. They also show that C-isotopes with JV=9~11, i.e., 15 C ~ 17 C, 
favor prolate deformation, which is natural as the new shell begins to be occupied after 
N=8. The 19 C nucleus with iV = 13 is shown most likely to favor oblate deformation. This 
suggests that the neutron number N=13 might locate in the latter half between the two 
closed shells and could be a manifestation of the new shell closure at N=16 instead of 
iV=20 as in the case of 0-isotopes||l|, It is thus desperately desired to have decisive 
experimental information on the signs of the deformations in heavier C isotopes. 

The neutron number dependence of the deformation in C isotopes was studied by using 



anti-symmetrized molecular dynamics (AMD) model in ref. [15]. They pointed out similar 
neutron number dependence for the deformation to the present results in neutron-rich C 
isotopes. Namely the neutron deformation changes from spherical in N=8, to prolate in 
iV=10 and then, to oblate in N=14, while the proton deformation stays always oblate 
independent to the neutron number. There is a difference in the case of N=16. The 
present deformed HF+BCS result shows a spherical minimum for 22 C, while the AMD 
model gives a triaxial shape for the neutron configuration. Since the present deformed 
HF+BCS model is performed assuming the axial symmetric deformation, the two results 
are not completely equivalent. It might be interesting to study further the deformation 
changes taking into account the tri-axial degree of freedom in the deformed HF+BCS 
model. 
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3 MAGNETIC AND QUADRUPOLE MOMENTS 
AND EFFECTIVE OPERATORS 



We next perform the shell model calculations for C isotopes with the effective interac- 
tions WBP10 in the (Op-Odls) configuration space |16| to study the magnetic and the Q- 
moments. The WBP10 interaction is designed to reproduce systematically the energy 
of ground state and excited states of stable sd shell nuclei. The energies and the spin- 
parities of the states near the ground states of odd C isotopes are tabulated in Table 2. 
It is interesting to see that two or three different spin states are almost degenerate in the 
odd C isotopes 15 C, 17 C and 19 C. These degeneracies are also expected from the results 
of deformed HF calculations, as we showed in the previous section. Several experimental 
efforts have been made to assign the spin-parities of odd C- isotopes 15 C, 17 C and 19 C. 
For this end, the magnetic moments and Q-moments will provide the most conclusive 
information. Calculated magnetic moments and Q-moments are given in Table 2. The 
effective spin g s — factor is taken to be 0.9g s (bare) for neutrons. This quenching factor is 
somewhat larger than the commonly adopted values 0.7—0.8 in stable nuclei. This differ- 
ence might be due to smaller effect of the second-order effects in the neutron-rich light 
nuclei |T7J]. In the ground state of 15 C, the calculated g factor is — 3.37/i/v which agrees 
well with the experimental one |g|=3.440±0.018 [|l8|l . The calculated values for the 3/2 + 
and 5/2 + states of 17 C are close to be — 0.514/ztv and —0.505^, while that of the l/2 + 
state is -2.82/i^. The empirical value |#(exp)|=0.5054±0.0025 p[ excludes the 1/2+ 
state from the ground state candidate, while the other two 3/2 + and 5/2 + states show 
good agreement within a few percent accuracy. There is a complemental experimental 
data of the selection rule on f3 decay from 17 C to 17 N which favors the spin 3/2 + as the 
ground state of 17 C. 

The effective charges e e ff(E2) are commonly adopted for the shell model calculations 
of Q-moments because of the limitation of the model space. In ref. pTJ, the polariza- 
tion charges e po i(E2) for the electric quadrupole moment are calculated by the harmonic 
vibration model and the isospin dependence is given by 

e pol (E2)/e = e eff (E2)/e - ~(1 - r z ) 
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= 3 x(r = 0) + - 2 x(r = D— 4- (-- X (r = 1) + ^(r = 0)j^— )r. 

(3) 

where x( r = 0) and \{r = 1) are the isoscalar (IS) and the isovector (IV) polarizability 
coefficients and ^P 1 is the ratio of TV and 75 components in the static nuclear potential. 
The polarizability coefficients are evaluated to be \( T = 0) =1.0 and x( T = 1) = —0.64 
by the harmonic vibration model and the ratio is taken to be -pg- = —2.6 from the 
empirical mean field potential strength. By substituting these values in Eq. (|3|), we obtain 

Z N — Z Z N — Z 

e pol (E2)/e = - - 0.32—^— + (0.32 - 0.65-— — )r z (4) 



In ref. ||21|| , a microscopic particle- vibration model was applied to calculate the polar- 
ization charges in C isotopes using HF and random phase approximations (RPA). This 
model gives state-dependent and isospin dependent polarizability coefficients. The aver- 
aged polarizability coefficients x( T = 0)=0.82 and \(t = 1) = —0.24 are found in 12 C to 
be smaller than those of the harmonic vibration model. The smaller \( r = 0) is mainly 
due to spreading of IS giant quadrupole resonances (GQR), while a large quenching on 
the x( T = 1) comes from a substantially small ratio of Skyrme interactions com- 
pared with the harmonic vibration model. The averaged polarization charges in active 
valence configurations are given by e po i(E2; n) / e = 0.53, 0.33 and 0.15 for neutrons and 
e po i(E2;p)/e =0.29, 0.16 and 0.05 for protons in 12 C, 16 C and 20 C, respectively. In 12 C, 
the calculated e po i values give e e //(n)=0.53 for neutrons and e e ff(p)=1.29 for protons, 
which agree well with the commonly used values e e ff(jp)— 1.3 and e e j/(n)=0.5 for light 
nuclei. The isospin dependence of the polarization charges might be papametrized from 
the values 12 C and 16 C to be 

,™w Z ,N-Z , jZN-Z, 
e pol (E2)/e = aj + b^— + (c + dj^—)r z (5) 

with 

a = 0.82, b = -0.25, c = 0.12, d = -0.25. (6) 

Eq. © gives e pol) (E2; n)/e=0.22 and e pol (E2;p)/e =0.07 for 20 C. A large difference 
between the calculated value and that from Eq. (0) is due to the effect of neutron skin 
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in 20 C. The small polarization charges in the very neutron rich nuclei are found to be 
important to explain the observed Q-moments of B-isotopes|2T|. 

We use the isospin dependent polarization charges in ref. [^TJ to calculate Q-moment of 
C isotopes in Table 2. In 17 C, the magnetic g— factors are essentially the same for the two 
configurations 3/2 + and 5/2 + and the calculated g-factors are close to the experimental 
value. The calculated Q-moments, however, are very different in the two configurations in 
magnitude and even in sign reflecting the different deformation of the two configurations. 
The neutron and the proton contributions for the Q-moment are \7 .\mb and 6Amb, 
respectively, in the 3/2 + state, while they are —5.5mb and — 3.8mb, respectively, in the 
5/2 + state. It should be noticed that the magnetic moment and Q-moment of the 3/2 + 
state in 17 C show large deviations from the single particle values ^(Schmidt) and, Q(s.p.) = 
— 37.7 e e ff(n)mb for 0d3/2 state, i.e., even the sign of these moments are different in the two 
calculations. The single particle Q-moment for the 0d^ 2 state is Q(s.p.)=— 53.9e e //(n)m6, 
which is much larger than the shell model prediction Q=— 9.3m6 ( Notice the standard 
value for e e //(n)=0.5 and the presently adopted value for 17 C is e e //(n)=0.33). However, 
the shell model Q-momemts are consistent with the prolate deformation for the 3/2 + state 
and the oblate deformation for the 5/2 + state which are suggested by the deformed HF 
calculations in Table 1. Thus the measurement of Q-moment will be the most decisive 
experiment to assign the spin and the parity of the ground state of 17 C and will provide 
experimental justification of the deformed HF+BCS and shell model predictions. 

The magnetic moments and the Q-moment of 19 C are given in Table 2. It is still 
under dispute whether the spin and the parity of the ground state of 19 C is l/2 + or 5/2 + . 
According to the shell model calculations, the lowest 3/2 + is also close to the lowest l/2 + 
and 5/2 + states in energy. The neutron and the proton contributions to the Q-moments 
are — 18.4m6 and —\A.7mb in the 3/2 + state, while they are — 0.6m6 for neutrons and 
l.Qmb for protons in the 5/2 + state of 19 C. The values for the 3/2 + state is consistent with 
the possible oblate deformation suggested in Table 1. On the other hand, the proton and 
neutron contributions have different signs in the 5/2 + state, and also very different from 
the single particle value for the 0d^ 2 state. This shell model results suggest the large 
configuration mixing in the lowest 5/2 + state of 19 C. It is seen clearly from Table 2 that 
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the magnetic moment and Q-moment are very different for each configuration in 19 C and 
will give decisive information on the spin assignment of the ground state of 19 C. 



4 GIANT AND PIGMY RESONANCES IN C ISO- 
TOPES 

The isovector (IV) giant dipole resonance (GDR) is the most well established collective 
mode throughout the mass table with large photoabsorption cross sections, exhausting 
most of the classical Thomas- Reich-Kuhn (TRK) sum rule (the energy weighted sum rule 
value) [03, ^3], 23|. As a microscopic model, we perform the shell model calculations for 



the dipole excitation mode in C-isotopes. The calculations take into account a model 
space of up to (l+3)hu excitations in 12 C, 13 C and 14 C including Os-Op-lsOd-lpOf shells. 
Other nuclei are studied in a model space of 1 Tiu excitation of Os-Op-lsOd-lpOf shells. 



The Warburton-Brown interaction WBP10jT6[ is used in this study with the model space 



(Os-Op-lsOd-lpOf). The center of mass spurious components in the wave functions are 
pushed up to higher excitation energies by adding a fictitious hamiltonian which acts 
only on the center of mass excitation |25|. In a restricted model space, there still remain 



some spurious components in the wave functions after the diagonalization of the model 
hamiltonian, especially when one uses HF or Woods-Saxon single-particle wave functions, 
instead of harmonic oscillator wave functions. In order to remove the effect of these 
spurious components on the transition strength, we use the effective transition operator 

A AT _ 7. y N N Z 

X = l = e Efe - -zj-YiY^h) = e- £ nY^h) - e- £ nY^h) (7) 

i i i 

in which the center-of-mass correction is subtracted from the IV dipole transition operator. 
The transition strength B(E1) is defined as 

B(El;uj n ) = J2\(n\d x = 1 \gs} | 2 (8) 
m 

where the matrix element is calculated between the ground state {\gs)) and the n-th 
excited 1~ shell model state (|n)) with the excitation energy huo n . In order to smooth out 
the discrete strength, the transition strength is averaged by a weight factor p(uj) as 



dB(El;u) 



■ ^B(El;u n )p(u - u n )<ko n (9) 



where 

1 172 

p(u-u n ) = — f trlW ( 10 ) 

The weight factor can be considered to simulate the escape and the spreading widths. The 
width parameter V is arbitrary taken as lMeV to draw a smooth curve of the transition 
strength. The oscillator length of the harmonic oscillator wave function is taken as b = 
{h,jmijj^) 1 l 2 =\. < o^ fm. It is known that the photoreaction cross section a is related with 
the transition strength B(El;u). The total photoreaction cross section a ini is written 
as|2§ 

r , 16tt 3 r E «™* dB(El;u) , , , 

*«* = J aduj= m^Jo u du duj 

The sum rule is a useful measure of the collectivity in GR. For the IV GDR, the energy 
weighted sum rule value is given by 

S(TRK) = |> n | <n|OfV) | 2 = ~~ = 14.9^ (MeV ■ fm 2 ) (12) 



neglecting the contributions of exchange terms. This sum rule (|i2D is known as the 
classical Thomas-Reiche-Kuhn (TRK) sum rule. The cross section <j int is then expressed 
as 

a int = ^S(TRK) = 60^ (MeVmb). (13) 

The cluster sum rule is referred to measure the adiabaticity between GR and Pigmy 
resonance PS]. Assuming the valence cluster with N2 and Z2 and the core with Ni and 
Zi, the cluster sum rule is given by 

s(cluster , = Jg.±(^ ; -W (14) 

v ; 2m4vr AA ± A 2 V ; 

where A\ = Ni + Z\ and A 2 = N 2 + Z 2 - 

The calculated results of averaged dipole strength (4) in C isotopes are shown in Figs. 
8 and 10. The non-energy weighted summed transition strength (NESR), the energy 
weighted summed transition strength (EWSR) and the total photoreaction cross sections 
Oi nt are tabulated in Table 3. The low energy strength below GDR region is compared 



with the cluster sum rule value fll4|) in Table 4. The photoreaction cross sections for C 
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13 C and 14 C are also shown in Fig. 9. The width parameter is taken to be T = 2 MeV in 
Fig. 9. 

In Fig. 8a, two calculated results of the transition strength in 12 C obtained within lhu) 
and (l+3)hu configuration space are shown. Although the peak energy of GDR appears 
at the same energy E x ~ 21 MeV, the integrated strength of the large space is 20% smaller 
than that of the small space because of the coupling to many-particle many-hole states. 
The experimental photoreaction cross section shows the GDR peak at E x = 22 MeVp4| 
that is close to the present calculation. The observed er^ value exhausts 64% of the TRK 
sum rule up to E x = 30 MeV ||24|| , while the calculated sums show the enhancement factor 
k for the sum rule, i.e., k=0.Q2 (0.29) in the lhu(3tkv) calculations. A large fraction of 
the strength is found to be in high excitation energy region up to E x = 140 MeV; k = 
0.62 up to 100 MeV and k = 0.86 up to 140 MeVjgg]. In ref. |7|], integrated cross section 
up to 35 MeV is found to be about 65% ~ 90% of the TRK sum rule, while calculated 
enhancement factor up to E x = 35 MeV is k = 0.01 in the (1+3) Tiuj space. In view 
of this, 25±10% of our calculated strength obtained within the (1+3) %uj space should 
be in higher energy region than E x = 35 MeV. We, thus, need to reduce the calculated 
cross section for 12 C by multiplying a factor 0.7 as shown in Fig. 9 to obtain reasonable 



agreement with the available experimental one in ref. [fPfl . The experiment data show the 
existence of a large fraction of the strength in higher energy region than E x =35MeV and 
suggest the importance of the coupling to many-particle many-hole states more than 3hu 
excitations. It is interesting to notice that the difference between the two calculations 
in Fig. 8a is only 200 keV for the peak energy, although the total cross section of 3hu> 
calculation is 20% smaller than that of lhu calculation. 

The calculated transition strength dB(El;uj)/duj and the photoreaction cross section 
u in 13 C ( 14 C) obtained by including up to 3hu excitations are shown in Figs. 8b and 
9b (Figs. 8c and 9c), respectively. We see appreciable cross sections in both 13 C and 14 C 
below 14MeV, while there is essentially no cross section in 12 C in the same energy region. 
The energy weighted sum (EWS) of the strength up to 14 MeV amounts to be 86% and 
66% of the cluster sum rule in 13 C and 14 C, respectively, as shown in Table 4. The GDR 
peak appears at 25~26 MeV in 13 C which is close to the experimental value at 24 MeV 
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p4j. The distribution of the observed photoreaction cross section is well reproduced by 



the present calculation within (1+3) %uo space. The integrated cross section calculated 
up to E x = 30 MeV amounts to be 98% of the TRK sum rule, while the observed one is 
71% of the TRK value About 30% of the calculated strength is in the higher energy 
region beyond E x = 30 MeV, which is similar to the case of 12 C. 

In case of 14 C, the GDR peak appears at E x = 28 MeV which is rather close to the ob- 
served one at 25.6 MeV[24J. Experimental mean energies are 18.3±0.4 MeV and 26.7±0.1 
MeV for T< = 1 and T>=2 states, respectively J24|, which are close to our calculated values; 
19.8 MeV for T< = 1 and 28.2 MeV for T>=2 states. The distribution of the photoreac- 
tion cross section is well reproduced by the present calculation. Observed photoreaction 
cross sections for T< and T> states summed up to 30 MeV are 88±12 MeV-mb and 37±8 
MeV-mb, respectively f24j], while the present calculation gives 119 MeV-mb for T< states 
and 94 MeV-mb for T> states. Experimental values are quenched compared with the 
calculated ones by factors 0.74±0.1 for T< and 0.40±0.08 for T> states. A factor 0.4 is 
multiplied for T>=2 states in Fig. 9. About 40±10% of the strength is found to be in 
the high energy region beyond E^. = 30 MeV. 

Calculated dipole strength for 15 C, 16 C, 17 C, 18 C and 19 C are shown in Fig. 10. Since 
the spin and the parity of the ground state in 19 C is not established yet experimentally, we 
calculate the dipole strength excited from two possible spin-parity states l/2 + and 3/2 + 
for the ground state. In case of 15 C, effects of skin are studied. The neutron lsi/2-orbit is 
obtained in a Woods-Saxon well to reproduce the experimental separation energy of 1.22 
MeV. The dipole strength enhanced about by 30% in the low energy region around E x 
= 5 MeV. The skin effect is rather moderate. The NESW and EWSR of the transition 
strength B(E1) and the total photoreaction cross sections are listed in Tables 3 and 4. 
The dipole strength below GDR region becomes substantial in these nuclei, i.e., the cross 
sections a int below hco=U MeV exhaust 7.8% for 15 C , 16.3% for 16 C , 13.1% for 17 C, 
11.6% for 18 C and 12.6% for l/2+ s . case of 19 C ( 14.1 % for 3/2+ s case of 19 C) of the TRK 
sum rule, respectively. These values correspond to 46.6% for 15 C , 81.3% for 16 C , 57.8% 
for 17 C, 46.4% for 18 C and 47.0% for 1/2+,. case of 19 C ( 52.4 % for 3/2+ s case of 19 C) of 
the cluster sum rule, respectively. The GDR peaks with the isospin T< are found always 



12 



at around E x = 17~19 MeV in these nuclei. On the other hand, the T> peaks appear more 
than 10 MeV higher in energy than the T< peaks, and smaller in peak height in heavier 
C-isotopes. In the extreme case of 19 C, the cross section a int of T> states becomes only 
13.5% of the TRK sum rule and there are only very small strength of T> states below 
E x =30MeV. 

There are peaks at rather low energies 10 - 20 MeV in 15 ~ 19 C. These energies are close 
to the unperturbed p— h excitation energy lftu;(H.O.)=41/A 1 / 3 MeV in the harmonic os- 
cillator model, but much lower than the systematic excitation energies of Giant resonances 
^(GDR^O/A 1 / 3 MeV. This low energy feature might be attributed to the effects of 
large deformations (see Table 1), which makes some unperturbed p-h 1~ states lower than 
those of spherical nuclei. The main GR part have two peaks in 15 C and 16 C. This can be 
also considered as the effects of the strong prolate deformation as is seen in Table I. In 
cases of 17 C, 18 C and 19 C, any clear two-peak structure is not seen in the strength distri- 
butions while the main peaks have large widths of T ~10MeV. The strength distributions 
show not much difference between the case of spin 1/2+ for the ground state (prolate 
deformation) and that of spin 3/2 + for the ground state (oblate deformation). As the 
strength around the peak up to 12-14 MeV exhausts about 50% of the cluster sum rule 
value in the heavier C isotopes, these regions may be interpreted as pigmy resonances. 
The strength is fragmented widely in the heavier isotopes , and the distinction between 
Giant resonance and pigmy resonance seems not very clear except for 15 C and 18 C. 

5 SUMMARY 

We have studied the ground state properties of C isotopes by the deformed Hartree-Fock 
+BCS model. The shallow deformation minima are found in the several isotopes. The 
prolate deformation is suggested to favor for 15 C and 17 C while the oblate deformation is 
most probable for 19 C. Both the prolate and the oblate minima appear in 16 C and 18 C to 
be almost degenarate in energy. 

The magnetic and the Q- moments of odd isotopes are investigated by the shell model 
calculations, and their configuration dependence is pointed out. It is crucial to obtain 
experimental information on the values of magnetic and quadrupole moments to determine 
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whether the deformation is prolate or oblate. In particular, it would be interesting to find 
out decisive information on the deformation of 19 C since this nucleus is a keystone to 
establish the new shell structure at N=16 in the C isotopes. 

We have also studied the Pigmy and GDR dipole strengths of C-isotopes by using 
shell model calculations in the large scale shell model (Op-lsOd-lpOf) space. We found 
that the excitation energies of GDR in 12 C, 13 C and 14 C show good agreement with the 
experimental data of the two isospin resonances, T> and T<. Moreover the calculated 
Pigmy strength below hu =14 MeV in 14 C is consistent with the experimental photore- 
action cross sections. In heavier C-isotopes than 14 C, the T< GDR has always a peak at 
around = 17~19 MeV, while the T> peak is more than 10 MeV higher in energy and 
much smaller in the cross section than the T< one. In these heavy isotopes, the Pigmy 
resonances are more pronounced than that of 14 C, having about 8~ 16% of the TRK sum 
rule values, which correspond to 50% of the cluster sum rule values . Future experimental 
effort is highly desirable to observe these Pigmy resonances to clarify the structure of drip 
line nuclei Q. 
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Table 1 The energy minima of the energy surface in the deformed HF calculations 
with the Skyrme interaction SkI4: (a) with the volume-type delta pairing interaction in 
Eq. (0), and (b) with the surface- type density- dependent pairing interaction in Eq. (§). 
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Table 2 Magnetic moments and quadrupole moments of C-isotopes. Shell model calcu- 
lations are performed with the Warburton-Brown WBP10 interaction. The effective spin 
g-factor for neutron is taken to be g s (ef f) / g s (bare)=0.9 in the shell model calculations. 
The effective charges for Q-moments are taken from the results of particle- vibration model 
based on HF+RPA calculations in ref. ||21|| . The experimental data of g— factors are taken 
from ref. [18 for 15 C and from ref. [19] for 17 C. The single particle Q-moment is calculated 
by using the harmonic oscillator wave function with the oscillator length b=1.64 fm. 
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Table 3 Non-energy weighted sum rule (NESR) and energy weighted sum rule (EWSR) 
values of El transitions in C-isotopes. Integrated photoreaction cross sections <Ji nt (MeV-mb) 
are also shown in the Table. Shell model calculations are performed with the Warburton- 
Brown WBP10 interaction. 
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Table 4 Low energy strength of electric dipole transitions in C isotopes. Energy- 
weighted sum rule (EWSR) values are compared with the cluster sum rule values a c i u . 
Shell model calculations are performed with the Warburton-Brown WBP10 interaction. 
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Figure 1: The energy surfaces for the ground state of even- mass C isotopes obtained by 
the HF+BCS calculations with a Skyrme interaction SkI4 together with the volume-type 
delta pairing interaction (upper panel) in Eq. ([[]) and the surface-type density-dependent 
pairing interaction (lower panel) in Eq. (0). 
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Figure 2: The energy surface for the ground state of even-mass C isotopes obtained by 
the HF+BCS calculations with a Skyrme interaction SkI4 together with the surface-type 
pairing interaction @ with the weak pairing strength V = —500 MeV-fm 3 for neutrons 
and —573 MeV-fm 3 for protons. 
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Figure 3: The energy surfaces for the ground state of even-mass C isotopes obtained by 
the HF+BCS calculations with Skyrme interaction SIII (upper panel) and SkM* (lower 
panel). The surface-type density dependent pairing interaction (§) is used for the BCS 
calculations. 
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Figure 4: One quasi-particle energy surface for 13 C on top of the BCS ground state of 
12 C. The SkI4 interaction with the volume-type delta pairing interaction (upper panel) 
and the surface-type density-dependent pairing interaction (lower panel) are used for the 
HF+BCS calculations. 
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Figure 5: One quasi-particle energy surface for 15 C on top of the BCS ground state of 
14 C. The SkI4 interaction with the volume-type delta pairing interaction (upper panel) 
and the surface-type density-dependent pairing interaction (lower panel) are used for the 
HF+BCS calculations. 
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Figure 6: One quasi-particle energy surface for 17 C on top of the BCS ground state of 
16 C. The SkI4 interaction with the volume-type delta pairing interaction (upper panel) 
and the surface-type density-dependent pairing interaction (lower panel) are used for the 
HF+BCS calculations. 
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Figure 7: One quasi-particle energy surface for 19 C on top of the BCS ground state of 
18 C. The SkI4 interaction with the volume-type delta pairing interaction (upper panel) 
and the surface-type density-dependent pairing interaction (lower panel) are used for the 
HF+BCS calculations. 
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Figure 8: Calculated B(E1) strength for 12 C, for 13 C and for 14 C with the use of the 
WBP10 interaction including up to 3hu excitations, (a) The solid curve shows the results 
with (l-\-3)huj excitations in 12 C, while the dashed curve gives those with 1 hu excitations 
only, (b) The solid curve includes both the results of T=l/2 and 3/2 states in 13 C, while 
the dashed curve gives the ressults of T=l/2 states only. The final states are J 7r =l/2 + 
and 3/2 + states, (c) The solid curve includes both the results of T=l and T=2 states, 
while the dashed curve gives the results of T=l states only. 



29 




Figure 9: Calculated photoreaction cross sectios for 12 C, 13 C and 14 C. The T> part of 
giant resonances is multiplied by quenching factors 0.7, 1.0 and 0.4 for 12 C, 13 C and 14 C, 
respectively. The shell model calculations are performed with the use of the WBP10 
interaction including up to 3hu excitations. 
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Figure 10: Calculated B(E1) strength for C-isotopes with N=9~13. with the use of the 
WBP10 interaction and the configurations of lhu excitations. The solid curve includes 
the results of both T< and T> states, while the dashed curve includes only those of T=< 
in all figures except in 15 C. (a) 15 C; The solid curve includes the results of both T< and 
T> states , while the dashed-dotted curve includes only those of T=< states . The dashed 
curve includes the effect of neutron skin for both T< and T> states, (b) 16 C. (c) 17 C. (d) 
18 C. (e) 19 C; The ground state is taken to be l/2 + . (f) 19 C; The ground state is taken to 
be 3/2+. 
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